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We derive a set of algebraic equations, the so-called multipartite separability eigenvalue equations. 
Based on their solutions, wc introduce a universal method for the construction of multipartite 
entanglement witnesses. Wo witness multipartite entanglement of 10'^ coupled quantum oscillators, 
by solving our basic equations analytically. This clearly demonstrates the feasibility of our method 
for studying ultra-high orders of multipartite entanglement in complex quantum systems. 
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I. INTRODUCTION 

Entanglement represents a fundamental quantum cor- 
relation between compound quantum systems. Since the 
early days of quantum physics this property has been 
used to illustrate the surprising consequences of the quan- 
tum description of nature [1, 2]. Moreover, entanglement 
plays a fundamental role in various applications and pro- 
tocols in quantum information science [3-5]. 

In multipartite systems a separable state is a statis- 
tical mixture of product states [6]. A quantum state 
is entangled, whenever it cannot be represented in this 
form. Various forms of multipartite entanglement are 
known [7 10]. The most prominent and non-equivalent 
forms of entangled multipartite quantum states are the 
GHZ-state [11] and the W-state [12], which have been 
generalized to so-called cluster and graph states [13, 14]. 
Another classification is given in terms of partial and full 
(or genuine) multipartite entanglement, for an introduc- 
tion see e.g. [4, 5]. Beyond finite dimensional systems, 
multipartite quantum entanglement in continuous vari- 
able systems turns out to be a cumbersome problem. 
Even in the case of Gaussian states, there exist multi- 
partite entangled states, which cannot be distilled [15]. 

The detection of entanglement is typically done via 
the construction of proper entanglement witnesses [IS- 
IS], being equivalent to the method of positive, but not 
completely positive maps. A witness is an observable, 
which is non-negative for separable states, and it can 
have a negative expectation value for entangled states. 
For different kinds of entanglement, different types of wit- 
nesses have been considered: bipartite witnesses [17, 19]; 
Schmidt number witnesses [20, 21]; and multipartite wit- 
nesses for partial and genuine entanglement [22-24]. A 
systematic approach for witnessing entanglement in com- 
plex quantum systems is missing yet. 

Recently, we considered the construction of bipartite 
entanglement witnesses with the so-called separability 
eigenvalue equations [19]. We have shown that the 
same equations need to be solved to obtain entanglement 
quasiprobabilities, which are non-positive distributions 



iff the corresponding quantum state is entangled [25]. 
Moreover, we have shown that the Schmidt number wit- 
nesses can be obtained by solving the related Schmidt 
number eigenvalue problem [21]. 

In the present contribution we derive a set of alge- 
braic equations, which yield the construction of arbitrary 
multipartite entanglement witnesses. For these so-called 
multipartite separability eigenvalue equations, we will 
study some fundamental properties, which uncover the 
structure of multipartite entanglement. Examples are 
given to witness partial and full entanglement in miilti- 
partite composed systems, for pure and mixed quantum 
states in discrete and continuous variable systems. 

II. MULTIPARTITE SEPARABILITY 
EIGENVALUE EQUATIONS 

In the following, we consider a composed Hilbcrt space 
y. = 'H\® ■ ■ ■ ®'Hn ■ It has been shown that, without loss 
of generality, we could assume that the individual sub- 
systems are finite dimensional ones [26]. Let us consider 
a partition Xi . . .Xk of the index set X = {1, . . . , iV}. A 
quantum state a is separable for the given partition, if it 
can be written as a classical mixture of product states [6] : 

0-= / AP{ai,...,aK)\ai,...,aK){ai,...,aK\, (1) 

with P being a classical probability distribution and 
Sxi:---:Xk being the set of pure and normalized separa- 
ble states. If a quantum state q cannot be written in the 
form of Eq. (1), it refers to as multipartite entangled. 

A multipartite entanglement witness for the given par- 
tition is a Hermitian operator W, with 

(TF) = ii{(jW) > 0, for all a separable, 

(1^) = tr(^W^) < 0, for at least one q. (2) 

Based on the Refs. [19, 27], it can be readily shown that 
any witness can be presented in the form 

W = fi,.......iAL)'^-U (3) 
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where L is a general Hermitian operator and the function 
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(^) denotes the maximally attainable expecta- 
tion value for separable states: 

fii:-:iK{L) = sup{(ai, . . .,aK\L\ai, uk)}- 

The suprcmum is taken over all |ai, . . . , ax) G Sij^-,...-,!^ . 

Hence, we can formulate a necessary and sufficient en- 
tanglement criterion being equivalent to the witness cri- 
terion: A quantum state g is entangled with respect to 
the partition Ii, . . .Ik, iff there exists a Hermitian op- 
erator L such that 

tr(^L)>/x,,.,x.(L). (4) 

This means that the mean value of L exceeds the bound- 
ary of niean values for separable states. A replacement 
L I— > —L leads to a similar entanglement criterion, but 
with the greatest lower bound (inf) instead of the least 
upper bound (sup): 

tr(^L) < inf{(ai, . . . , aK\L\ai, ax)}- (5) 



For both entanglement criteria, we have to solve the 
following optimization problem for an observable L: 

G = (ai, . . . , aK\L\ai, . . . , uk) — > optimum, 

C = {ai,...,aK\ai,...,aK) -1 = 0, (6) 

where G represents function to be optimized, and G is 
the normalization condition. For such an optimization 
problem, we can apply the method of Lagrangian multi- 
pliers. In our case, the optimization conditions is 



d{aj\ d{aj\ 

where g is the Lagrangian multiplier and is the null 

vector in the subspace given by the partition Xj. The 
partial derivatives of G can be computed as: 



dG _ g(ai,...,ajv|L|ai,...,ajv) _ ^ 
d{aj\ d{aj\ 



trxi ...trx^ (^L\ai, . . . ,aK){ai,. . . ,aK\j 
d{aj\ 



=trxi . . .trx^._itrxj.,.i . . .trx^^ (j^ [\ai, . . . ,aj_i)(ai, . . . ,aj_i| (g) Ixj (8 \aj+i, ■ ■ .,aK){aj+i, . . . ,aif |]) \aj) 



-Ijai,...,aj-i,aj+i,...,aK I'^j) ■ 



(8) 



The case L = 1 yields the derivatives of C. Let us also 
note that we assiimed that the indices of the sets X, are 
ordered in a form that all elements of Ij are larger then 
the elements of Iji for j > j'. This assumption is justified 
by the fact that one can employ, without loss of gener- 
ality, a permutation of the Hilbert spaces "Hi . . .T-Ln to 
order them in the required form. 

The Euler-Lagrangian optimization condition in 
Eq. (7) can be reformulated for all j = 1, . . . , X as 

= iai,...,a,_i,a, + i,...,a,c|aj) ^ gWj) , (9) 

where all eigenstates are normalized ones, {aj\aj) = 1. 
In addition, we may evaluate the value of g. We can do 
this by multiplying Eq. (9) with {aj \ . This yields 

9 —{<^j\-^ai,...,aj-i,aj+i,...,aK I'^i) 

=(ai, . . . ,aK\L\ai, ■ ■ ■ , o-k) = Goptimum- (10) 

Hence, the Lagrangian multiplier corresponds to an op- 
timal expectation value of L for separable states. In con- 
clusion of this derivation, we get an algebraic problem 
whose solutions give all optimal expectation values. 
Definition (MSEvalue equations). The equations 

La,,...,aj-i,aj+i,...,aK\aj) =51%) ioT j = 1,...,K 



are defined as the first form of the multipartite separa- 
bility eigenvalue (MSEvalue) equations. The value g is 
denoted as the MSEvalue of L, and the product vector 
|ai, . . . , gk) is the corresponding multipartite separabil- 
ity eigenvector (MSEvector). 

As a final conclusion from this derivation we get 

fI^■.■■■■.IK{L) = sup [g: g is MSEvalue of i} , (11) 

and condition (5) is given by the infimum of all MSEval- 
ues. This means that all multipartite entanglement wit- 
ness can be constructed from the solutions of the MSE- 
value equations, 

W = sup{g}i - L (12) 

In case we consider finite Hilbert spaces (dim'H < oo), 
we can replace sup and inf by max and min, respectively. 

The derived MSEvalue equations play a fundamental 
role for multipartite entanglement tests. They give the 
possibility to construct arbitrary entanglement witnesses 
on the basis of the solutioii of an algebraic eigenvalue 
problem of an observable L. Numerical and analyti- 
cal methods - originally developed to solve eigenvalue 
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problems - can be applied to handle the multipartite en- 
tanglement problem in quantum physics in a systematic 
way. Before we apply our method, let us formulate some 
fundamental properties of the MSEvalue equations. The 
proofs are given in Sees. A-C of the Appendix. 

Proposition (Second form of MSEvalue equations). 
The Hermitian operator L has the MSEvalue g for the 
MSEvector |ai, . . . jUk), iff it fulfills the second form of 
the MSEvalue equations 

L\ai, ...,aK) = g\ai, . . . , ox) + |x), 



,aK\x) 



for all I a;) e 

□ 



with (ai, . . . , a^-i, a;, aj-i-i, . . 

and j = l,...,K. 

This proposition transforms the coupled system of 
eigenvalue equations, which has been defined in the first 
form of the MSEvalue equations, into a single, but per- 
turbed eigenvalue problem. This second form yields sev- 
eral implications. For example, if an eigenvector is a 
product vector |ai, . . . , a^), it is also an MSEvector with 
Ix) = 0. In addition, we also conclude that the operator 
L yields a true entanglement witness, cf. Eq. (2), ifi^ the 
eigenspace of the largest eigenvalue, does not contain a 
product vector. 

Proposition (Transformation properties). A Hermi- 
tian operator L has a MSEvalue g for the MSEvector 
|ai, . . . , a^)- Then the operator 



AiH-A2£ 



(t/i O • • • Uk) 



with Ai,A2 S M \ {0} and C/j being unitary opera- 
tions acting locally on the partition Ij, has the MSE- 
value g' = Xi + \29 and the MSEvectors \a'{, . . . ,a'^) = 
Ui® ■ ■ ■ ®Uk |ffli, • • • ,0^)- □ 
This transformation allows us to consider a whole class 
of witnesses, by solving the MSEvalue equation for a par- 
ticular operator L. In addition, the shifting of L to L' al- 
lows us to consider positive semi-definite operators only. 
Note that the invariancc of the MSEvahies under local 
unitaries is of particular interest for the quantification of 
multipartite entanglement, see e.g. [4, 5]. 

Proposition (Cascaded structure). The non-zero so- 
lutions of an TV -I- 1-partite operator L' = |V-')(V'I si's 
identical to the solutions of an iV-partite operator L = 

tVN+lL'. □ 

This property is quite surprising. It shows us that 
all possible entanglement witnesses - based on positive 
semi-definite operators L - of an A/'-partite system can 
be constructed by a few simple entanglement witnesses in 
a A -|- 1-partite system, L' = An arbitrary rank 

of L can be achieved by choosing a state with the 
same Schmidt rank for the bipartition Ii = {1, . . . , N} 
and I2 = {A -|- 1}, cf. Sec. C in the Appendix. 



III. EXAMPLES 

In the following we apply our method to analytically 
derive multipartite entanglement tests. First we may 
consider witnesses for prominent examples of states in 
a three qubit systems. In a second step, we apply our 
method to get a multipartite entanglement test in a com- 
plex continuous variable system. 

Let us consider a generalized tripartite W-state 

IV'w) = Ai|l,0,0) + A2|0,l,0)+A3|0,0,l), (13) 

with |Aip-|-|A2p-|-|A3|^ = 1, which defines the observable 
L = \tpw){'>Pw\- In Sec. D of the Appendix, we solve the 
MSEvalue equation of L. This gives 

/{i}:{2,3}(i) =max{|Ai|MA2|' + lAsI^}, 
/{2Hi,3}(i) =max{|A2|MAi|2 + |A3|2}, 
f{3}:{i,2}{L) =max{|A3|^,|Ai|^ -|- IA2P}, 
f{i}:{2}:{3}{L) = max{| Ai |M A2 |M A3 5o}, (14) 



with 



90 = 



|AinA2p|A3p 



|Ai|2 + |A2|2 + |A3|2)2 - 2(|Ai|4 + {X^]^ + jAal*) ' 



Hence, we can formulate the following multipartite en- 
tanglement conditions: A quantum state p is partially 
entangled, if (V'wIpIV'w) > /{i_}:{2}:{3}(-^)- The corre- 
sponding entanglement witness is 

Wpart = max{go, |Aip, IA2P, lAgp}! - |Vw)(V'w|. (15) 

A quantum state p is fully entangled, if (V'wIpIV'w) > 

niax{/|i};|2,3}(-^)'/{2}:{i,3}(i'),/{3}:{i,2}(^)}- The cor- 
responding entanglement witness is 

Wfuu = max{|Ai|' + |A,f : i ^ j}i - |^w)(V'w|. (16) 

In Fig. 1, we apply the considered witness to study the 
entanglement of a noisy W-statc. 

In a second step, a generalized GHZ-state is given. 



IV'GHz) =ko|0,0,0)-Fki|1,1,1), 



(17) 



together with |«;o|^-|-|ki|^ = 1, which yields an observable 

L = |'(/'GHz)(V'GHz|- From Sec. E of the Appendix, we get 
the maximal MSEvalues: 

/{l}:{2,3}(-^) = /{2}:{l,3}(-^^) = /{3}:{1,2}(-C') 

=/{i}:{2}:{3}(^) = max{|KoP, (18) 

Hence, a state p is genuinely tripartite entangled, if 
(V'GHz|/5|'0GHz) > max{|Ko|MKiP}, see Fig. 2. Note 
that the corresponding witness 

W = max{|Ko|', - |V'GHz)(V'GHz| (19) 
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FIG. 1: (Color online) The entanglement test in Eq. (4) for a 
W-state mixed with white noise, p = p|l + (1 — p)l^w)(i/'w| 
Ai = A2 = A3 = is plotted for < p < 1. The boundary 
for partial or full separability is | or |, respectively. The 
expectation value (L) exceeds the boundary for full or partial 
separability as long as the mixing parameter is p < |§ or 
P < respectively. 
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FIG. 2: (Color online) The expectation value (L) for a GHZ- 
state mixed with white noise, p = + (1 — p)|i/'ghz){V'ghz| 
(ko = Ki = ^), is plotted. The lower boundary for par- 
tial and full entanglement is |. The expectation value (L) 
exceeds this boundary - which automatically implies genuine 
entanglement - as long as the mixing parameter is p < j^. 



cannot discriminate between partially and fully entangled 
states, see Eq. (18), which is possible for the generalized 
W-state projection in the previous example. 

As a proof of principle, we are going to test multipar- 
tite entanglement of a continuous variable system. Our 
considered example is a system of N coupled harmonic 
oscillators. The observable we are using to verify entan- 
glement is the total energy of this system, L = H, 

where 7 denotes the coupling strength of the interaction, 
fj the position and the momentum operator. Let us 
note that we considered an even more general case in 
Sec. F of the Appendix. For the partition Ii, . . . ^'I■K■, we 
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FIG. 3: (Color onhne) The energy difference A = ((H) - 
i5[I])/ (fftw) between quantum states and the ground state 
of the system is plotted depending on the interaction strength. 
The boundary Aiooo corresponds to the minimal attainable 
energy for fully separable states. All quantum states having 
an energy below A5 (A''i = • ■ • = = 200) reveal entangle- 
ment between more than 200 oscillators. The quantum states 
in the range below A2 (A'^i = N2 = 500) exhibit entanglement 
distributed over more than 500 subsystems. 

get the smallest MSEvalue of the Hamiltonian as 

+ f^W^^^.)^ (21) 

where Nj = \Tj \ is the number of subsystems in Tj . The 
resulting witness reads as 

Wi,,...^i^=E[Ii,...aK]l-H- (22) 

In the special case K — 1 {Ii — I) we get the true 
ground state energy of the system, 

E[I]^lf^[f7^,{N^l)^l). (23) 

In case of full separability, K — N {Ij — {j}), we have a 
minimal energy of: 

E[{1}, . . . , {N}] = l + [N-l]^. (24) 

In Fig. 3 we plotted the corresponding entanglement test 
based on Eq. (5), for TV — IQ^ interacting oscillators. 
For the witnessing by the total energy iJ, no information 
about the structure of the quantum states is needed. 

IV. CONCLUSIONS 

In conclusion we have derived an algebraic set of equa- 
tions to construct arbitrary entanglement witnesses. We 
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studied some fundamental properties of these equations. 
For example, they are invariant under local unitary trans- 
formations and have a cascaded structure. The latter al- 
lows to deduce all entanglement witnesses from elemen- 
tary projections. Our method enables us to use all known 
proc;cdurcs for solving eigenvalue problems to construct 
entanglement witnesses. We applied our method to ana- 
lytically identify full and partial entanglement of general- 
ized, noisy GHZ- and W-states. Moreover, we witnessed 
multipartite entanglement for a system of 10"^ interact- 
ing oscillators, by analytical computation of the energetic 
boundaries of separable states. This demonstrates the 



feasibility of our method for studying ultra-high orders 
of multipartite entanglement, which are of fundamental 
interest for understanding the transition from the micro- 
scopic to macroscopic world. 
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Appendix A: Proof — Second Form of MSEvalue equations 

Let us assume that g and |ai, . . . , ax) are a MSEvalue and a MSEvector of L, respectively. This means they fulfill 
-£ai,...,aj_i,aj+i,OA: 1%) = for all J = 1, ... , K. The result of the mapping of L acting on a vector |ai, . . . , ax) can 
be decomposed in one part parallel to the input state and another part |x) perpendicular to the input state 



L\ai,...,aK) =g'\ai,...,aK) + \x) 



(AI) 



Note that g = (ai, . . . , aK\L\a\, . . . , qk) = g' ■ Let us consider an arbitrary |a;) € 0ieij. = !> • • • i We get 
the following projection for the map L — gl: 

0=g-g= {X\{L - gi)au-,aj-ua, + uaK\aj) 

={ai, . . .,aj-i,x,aj+i,aK\ L\ai, . . . ,ax) -g{ai,.. . ,aj-i,x,aj+i,aK\ai, . . .,aK) 

^ y 

=g\ai,...,aK) + \x) 

=g{ai, aj-i,x, aj+i, |ai, . . . , gk) + (ai, . . . , aj-i, x, aj+i, aK\x) 

- g{ai, a^-i, x, a^+i, a^loi, . . . , Uk) 
=(ai, . . . ,aj_i,x,aj+i,aK\x)- (A2) 
Hence, the MSEvalue equations imply that for the solution holds the form 

L\ai, ....qk) =g\ai,. ..,aK) + \x), (A3) 

where |x) has the property given in Eq. (A2). The other way around, we get from these equivalent transformations, 
that Eq. (A3) implies that g and |ai, . . . , uk) are a MSEvalue and a MSEvector of L, respectively. 



Appendix B: Proof — Transformation Properties 

Let us assume that g and \a\,. . . ,ax) are a MSEvalue and a MSEvector of L, respectively. Hence, Eq. (A3) is 
fulfilled. Now, we consider an operator L' (Ai, A2 e M \ {0}): 



- \1 
Uk) 



X1I + X2L 



•Uk) 



(Bl) 



where Uj are unitary transformations acting on (^^gx, (j = 1) -?'')• We get for the transformed input vector 
L'\a[, a'x) =L' (ul\ai) ^---^ U^ai)) 

= (^C/i (g) • • • O Uk^ (Ai|ai, . . .,aK) + X2g\a.i, . . . ,aif) + \x)) 

=(Ai + X2g)\a[, . . . ,a'^) + (j/i • • • Uk)^ \x) = (Ai + A25)|ai, • • • ,<) + Ix'), (B2) 
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where \x') obviously fulfills Eq. (A2) for \a[, ... , a'^) and arbitrary |x') = Uj\x) G Hence, Eq. (B2) is the 

second form of the MSEvalue equation for L' with the transformed solutions: 

= Ai + A25 and \a[, . . . ,0'^^) = f/||ai) f/j^lax). (B3) 

Let us note that an inverse transformation of L' yields the inverse implication, 

L={u^(^---(^Uk) -yl + y^V (j7i0---0C/k)\ (B4) 

Appendix C: Proof — Cascaded Structure 

Wc consider a rank one operator L' = Let us consider the Schmidt decomposition [3] of l'^), with respect 

to a bipartite decomposition {!,..., A''} and {N + 1}, as 

r 

\i') = Y.^j\i'j)®\^j), (ci) 

with Xj e R>o (Schmidt coefficients), r (Schmidt rank), orthonormal € ^f^^iHi, and orthonormal \zj) € Hn+i- 
The MSEvalue equation in the first form for the N + 1-th subsystem reads as 

r 

g\aN+i} =^'al,...,a^,kJV+l) = {i)\ai, . . . , gn+i}^ Xj {ai, . . . , aN\i^j)\zj) , (C2) 

with the abbriviation for the MSEvalue g = I7P (7 = {ip\ai, . . . ,aj^^i)). Here, we are only interested in solutions 
with 5^0. Prom Eq. (C2), we get that the solution |ajv+i) must have the form 

\aN+i) = 2^ — \Zj). (C3) 

We may insert this result into the MSEvalue equation for the i-th party {i = 1, . . . , N): 

gWi) =^ai,...,a._i,a.+i,...a„ + i l«»> 



=(V'|ai, . . ■ ,ajv+i) y^^Xj{aN+i\zj) ((ai, . . . ,ai_i| ig) • (g) (aj+i, . . . ,ajv|) \i>j} 
^ X^ 

=7^ — ((ai,. . . ,ai_i| (g) • (g) (oi+i, ...,aN\) \ipj){ipj\ai, . . . ,ajv) 



.•=1 ^ 

= ((ai,...,ai_i| (g) • (g (ai+i,...,ajv|) A||V'j)(V'jl j (|ai, . . . , ai_i) ig) la^) ig) |oi+i, . . . ,ajv)) 
where we used the abbreviation 

r r 

L = tr^+i|^)(V'| = ^ A,-A,.tr;v+i(k,-)(^/|)|V,-)(V',"l = E ^'l^^'^^^'l- (^5) 
i,i'=i j=i 

Let us note that L is a positive semi-definite operator which has a rank r. 

Appendix D: W-state witness 

Now, we consider the operator L = |V'w)(V'w|i where |'0w) = Ai|l,0,0) + A2IO, 1,0) + AajO, 0, 1) is a generalized 
W-state. The reduced operator - with respect to the third subsystem - is 

L'=tr3L= |^')(V''l + |A3nO,0)(0,0|, with |V') = Ai|l, 0) + A2IO, 1). (Dl) 
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We may start with the calculation of the maximal MSEvaluc for partial separability in the case {3} : {1,2}. Since 
and |0, 0) are perpendicular vectors in "Hi ® 7^2, we immediately get 

5max = max{|A3|2, (^'IV') = |Ai|2 + jAal^}. (D2) 

Similarly, this yields the maximal MSEvalues for the other decompositions: 5max = inax{|Ai|^, |A2p + lAa]^} for 
{1} : {2,3} and ^^ax = max{|A2|2, jAip + [Aa^} for {2} : {1, 3}. 
For the full separable case {1} : {2} : {3}, we can easily get the obvious solutions of L': 

(7= IA3P for |ai,a2) = |0, 0), g = |Aip for |ai,a2) = |1,0), (D3) 
g = IA2P for |ai,a2) = |0, 1), and g = for |ai,a2) = |1, 1). 

A convenient parametrization of all other possible MSEvcctors is 

M = « for a„a2 . C \ {0}. (D4) 

\/l + |air Vl + |a2r 

Here, we use the MSEvalue equations in the second form L'|ai, 02) = fflai, ci2)+|x)- Since we require that (a;i,a2|x) = 
and (ai,X2|x) = for all \xi) {i = 1,2), we have 

\/i + |air vi + F2r 

This leads to the MSEvaluc equation in the second form as 

|A3||0,0) + A2(A;a2 + Atai)|0, 1) + Ai(A;a2 + A*ai)|l, 0) + 0|1, 1) 
=(<? + halal)\Q,0) + (<?a2 - hal)% 1) + {gai - hal)\l,Q) + (501^2 + h)\l, 1), (D6) 

where we multiplied the whole equation with the normalization constant -y/l + |aip-\/l + |a2p. We may decompose 
ai = rjc"^' for i = 1, 2 in polar coordinates. The individual basis components can be rewritten as 

lAal" =5 + /iVir2 (D7) 

A^(A^*r2 + A'i*ri) =<?r2 - h'n (D8) 

A'l (A'2*r2 + Ai*ri ) =gn - h'r^ (D9) 

=grir2 + h', (DIO) 

where we introduced h' = /ie~'('^i+*'=' and A- = AjC"'^' {i = 1,2). Since g = (ai, a2|I/'|ai, a2) € M, it follows from 
Eq. (DIO), that h' € K. Eqs. (D8) and (D9) can be combined to 



A'iA'2* + h' |A'i|^ -g J \ rij ~ \ 



(Dll) 



Note that this equation has a real solution (ri,r2), iff lm{X'2X'*) = or, equivalently, Xi^X'* = X[X'2 ■ More generally, 
this equation has a solution, iff 

= det f l^'^X 7 1 t = i\Xif - 9){\X2\' -g)- |A'2Ai* + hf. (D12) 



,X[X'2*+h' \X[\-'-g 
In addition, we may resolve Eqs. (D7) and (DIO) as 



= r,r2 = ^^h'^= g{g - \X,n (D13) 

Hence, we get quadratic equations (D12) and (D13) in g and h' . The difference of both equations gives the relation 

between g and h' as 
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This relation, we may insert into Eq. (D13) which results in 



= 5 

Omitting g' = 0, we get the root of this equation as 



(|A3|^-|Ai|^-|A2|^) 
4|Ai|2|A2r 



2\2 

2 



Aar • (D15) 



4|Ai|^|Ao|-'|A;i|^ 

ffo 



2(|Ai|2|A2|2 + |Ai|2|A3|2 + lAsPIAsP) - |Ai|4 - |A2|4 - lAsI^ 

4|Ai|2|A2|2|A3|2 



(D16) 



(|AiP + |A2|2 + |A3|2)2 - 2(|Ai|4 + |A2|4 + jAal^)' 
Finally, we conclude for the case of full separability that the maximal MSEvalue gmax of L is 

5^a. = max{|Ai|MA2|MA3|2,ffo}, (D17) 

where we have taken the solutions in Eq. (D3) into account. It is of importance to mention that from the calculation 
of the MSEvector of Qq follows, that we have to fulfill the requirements |Ai|^ + IA2P > IA3P, IA3P + IA2P > |Aip, and 
|Aip + |A3|2> IA2P. 

Appendix E: GHZ-state witness 

Let us consider the operator L = |V'GHz)(^i'GHz|, with I'^ghz) = «o|0, 0, 0)+k;i|1, 1, 1) being a generalized GHZ-state. 
The reduced operator is 

L' = tr3L= |«or|0,0)(0,0| + |Kini,l)(l,l|. (El) 

First, wc may consider a decomposition {3} : {1, 2}. Since L' is already given in a spectral decomposition, we get the 
non-zero MSEvalues |k;oP and |ACip. Similarly, we can argue for the partial separability with respect to {2} : {1,3} 
and {1} : {2,3}. We get that the maximal MSEvalue, ^max, for partial separability is 

5max = max{|«;o|^, Iki^}- (E2) 

Second, wc consider full separability {1} : {2} : {3}. Due to the reduction from L to L' , we have to solve the 
MSEvalue equations of L' as given in Eq. (El). It is easy to check that is an MSEvector for the MSEvalue 
{i = 0,1). For all other MSEvectors, we use the same parametrization as given in Eq. (D4). Now, the MSEvalue 
equation read as 

fix I , ^ |Kop|0) + |Ki|'|ai|^a2|l) |0)+a2|l) 

^ai|a2)=g|a2)^ I2^ /i^i 12 = ^ /1 ^ 1 12 ^^^^ 

(l + |ai| )vl + |a2r \/l + |a2r 

f U \ .1. \ ^ |AConO) + |^ip|a2|V|l) |0)+ai|l) 

K2m) = 9\ai) ^ —— — _^ I |2 ~ =g /rxi= if ^ ' 

These equations are fulfilled, if |ai| = \a2\ = \ko\/\ki\ for the MSEvalue g = |koP|kiP/(|«oP + l/^iP)) see the 

components in |0) direction. Due to the fact that this g is smaller or equal to we get the maximal MSEvalue 
for full separability as 

flmax = max{|KoP, |ki|^}- (E5) 

Appendix F: Energy Witnesses 

Let us consider the Hamilton operator H oi a. system of N harmonic oscillators. 

H = -^V^M-W+lx'^Gx, (Fl) 
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where wc used the position operator vector x = {xn)n=i ^^'^ momentum operator vector p = 5V = 5 {dx„)^-i- The 
masses of the individual oscihators are given by 



M = diag(mi, . . . ,mAr), 



(F2) 



and the total potential energy reads as V{x) = x'^Gx, with: G G M.^^'^ (real valued); G = G^ (symmetric); and 
G > (positive definite) . The well-known ground state is given by the normalized, Gaussian wave-function 



detV 



1 



x'^Vx 



with V = T' 



In order to get V, we need the first and second derivative: 

\/^po{x) ^ {~Vx)i;o{x) and \/ (»\/ipo{x) = {-V + Vx (»Vx)tpo{x) 
The resulting eigenvalue problem is 

"I r fe2 1-1 



-Tr{M-'V) 



xipo{x) = Eoipo{x), 



(F3) 



(F4) 



(F5) 



where the here used trace operation, Tr, is acting on M"'^^"'^ matrices. The eigenvalue equation is solved by the ground 
state tlJo{x) for 







The minimal eigenvalue is 



2 2 hi \ 



^o = ^Tr 



(F6) 



(F7) 



As an intermediate step, it is useful to show that a minimal energy in our case is always obtained for states with 
zero mean values for x and p. We have an energy given by 



{H) =tJm-\-%WW^)] + It.{G{xx^)) 



=TrQM-i [{{Ap){Apf) + {p){pf]^ + ^Tr(G [{{Ax){Ax)'^) + {x){x)^]) 
--l{pfM-'(p) + l{xfG{x) + jTr[M-iCovar(p,pT)] + lTr[GCovar(x,fT)] . 



(F8) 
(F9) 
(FIO) 



A local displacement operation gives the smallest value for the case (p) = = {x), when considering the same 
covariance matrices. 

Now, we may study an arbitrary decomposition of X = {1, . . . , N}: 



K 



(I,)f=i: I,- CI, I,-^0, I,nX,. =0(fori^/). IJ = ^• 



(Fll) 



Wc define the position and momentum operator for the corresponding subspace as X(j-) = {xn)neij ^^^d Pij) = 
j{9x^)neijj respectively. The block-matrices M(^j_ji) and G(^jji) are similarily obtained from M and G, by ignoring 
all rows or columns which are not in Xj or X,'. Note that Mf^j^ji) = for j ^ /. A product wave function of this 
decomposition (2j )jLi reads in position representation as 



K 



Hence, we may write the energy operator as 

K 

2 



1 ^ _ 1 ^ 



(F12) 



(F13) 
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Now the MSEvalue equations read for a particular choice of J = 1, ... ,K as 



■0(1) >•■■,'/'( J-l),'/'(J+l)v>'/'(iV) 



(F14) 



together with the partially reduced operator 



V'(i)i---.V'(j-i).V'(j+i).--->V'(N) 



(J) 



3^J 



+ 9^(J)G(J,J)X(j) + 



(J), 



(F15) 



jjtjj'jtj 



where we used (p) = = {x) . Now the MSEvalue equation for the partition Ij is an ordinary eigenvalue problem of 
a quadratic (reduced) Hamiltonian: 



h"^ _ 1 

-^^Jj)^{Ij)^{J) + 2'^Jj)G{J,J)X(j) 



V'(J) (a^(J)) 



(F16) 



Prom our initial solution of the ground state of the full Hamiltonian, we obtain the ground state of the reduced 
Hamiltonian being the MSEvector in position representation, 



V'(J)(a^(j)) 



detVj 



exp 



Jj)VjXiJ) 



, forV, = iM;0) 



M 



1/2 
(J, J)' 



with the cardinality Nj oil j, Nj = \Tj\. This yields the corresponding minimal MSEvalue as 

K K ^ 



Eo[(2:.-)f=i]=E^o[^i] = E> 



i=i 



i=i 



Let us note, that we get for K = 1 Eq [{I)] = Eq, cf. Eq. (F7). 
The example in the paper has the following properties: 

M =mdiag(li, . . . , In), 



G = [niLj^ + 7V7] diag(li, . . . , 1a,) - [Nj] nj 



«=^(l^)^=l 



N 



(F17) 



(F18) 



(F19) 
(F20) 

(F21) 
(F22) 



Note that each component is a 3 dimensional matrix itself and the potential energy is y(ri, . . . , fjv) 



^ [(fj)^i] G [(fj)^i] . For a given Ij, we have 



(1.). 



n 



-1/2- 



2v^ 

h 

2\prn 



TrM [mw2 + [diag(li, . . . , l^v) 



^0) "0) 



+ [mw2^7(iV- |J,-|)] na)n5)) 



Tr \mu? + ]<S^\ 



1 + 



7iV 



[3|X,|-3] + 



diag(li)iei,- - %) ^^•) 

7(iV-|I,|) 



1 + 



+ [ma;2+7(7V-|Z,.|)]^fi(,)n^.^) 



(F23) 



(F24) 
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where the second line of Eo[Ij] represents the spectral deco mposition of Gqj) . 
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